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Abstract 

We write down a generator program for the hadronic production of the double-heavy baryons 
H cc , Hfe c and according to relevant publications. We name it as GENXICC and we test it by 
comparing its numerical results with those in references. It is written in a PYTHIA-compatible 
format and it can be easily implemented into PYTHIA. GENXICC is also written in modulariza- 
tion manner, with make, a GNU C compiler, one may apply the generator to various situations 
or experimental environments very conveniently. 
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Program summary 



Title of program: GENXICC 
Version: 1.0 (January, 2007) 



Program obtained from: http://www.itp.ac.cn/~zhangzx; CPC Program Library 



Computer. Any LINUX based on PC with FORTRAN 77 or FORTRAN90 and GNU C 
compiler being installed 

Operating systems: LINUX 

Programming language used: FORTRAN 77/90 

Memory required to execute with typical data: About 2.0 MB. 

Distribution format: tar gzip file 

Nature of physical problem: Hadronic production of a double-heavy baryons: H cc , E bc and E bb . 

Method of solution: The production of the double-heavy baryons is realized by producing 
a binding double-heavy diquark either (QQ')[ 3 Si}3 i6 (Q,Q' = b,c) or (QQ') P So] 3,6; which 
is in color anti-triplet 3 or color sextuplet 6 and in S-wave triplet or singlet configuration 
respectively, and then by absorbing a proper light quark non-perturbatively. For the 
production of the various double-heavy baryons S cc , H^ 1 and 5 66 , the 'gluon-gluon fusion' 
mechanism, being the most important, is written precisely in the generator, but two 



1 In fact, there are two kinds of states for S(, c , i.e., one is that the inside b and c are symmetric in 'flavor 
space' and the other is that b and c are anti-symmetric in 'flavor space' similar to the case for the baryons 
A and S . Let us call them as Si& c for symmetric one and "E.2bc for antisymmetric one when we need to 
distinguish them. Due to the electromagnetic interaction between the quarks for instance, the two kinds 
of states may have different masses (degeneracy broken) . 
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additional mechanisms, i.e. the 'gluon-charm collision' and the 'charm-charm collision' 
ones, only for S cc (S+ or 2+ + ) are written. Furthermore, all the mechanisms are treated 
consistently within the general-mass flavor-number (GM-VFN) scheme. Specially, to deal 
with the amplitude and in order to save CPU time as much as possible, the 'improved 
helicity-approach' is applied for the most complicated gluon-gluon fusion mechanism. The 
code with a proper option can generate weighted and un-weighted events accordingly as 
user's wish. Moreover, an interface to PYTHIA is provided to meet ones' needs to generate 
the 'complete events' of Hqq', i.e. to do the 'showers' of the partons appearing in the initial 
and final states of the subprocess, and the hadronization for final obtained 'showers' etc. 

Restrictions on the complexity of the problem: In GENXICC, the approach to the hadronic 
production in terms of a 'complete calculation' via the production of a binding diquark 
state either (QQ)[ 3 Si]^ or (QQ)[ 1 S }6 (Q = c,b) for S cc and E bb production, and via that 
of a binding diquark state of (bc)[ 3 Si}^ or (bc)[ 1 S \^ or (bc)[ 3 Si} 6 or (6c)[ 1 S , ]6 for E bc is 
available, but the contributions from the other higher Fock states of the diquark states are 
not involved. Considering the needs of comparisons and applications in most cases, three 
mechanisms and their consistent summation for the hadronic production of S cc are available. 
But for most purposes and applications to the baryons S&& and S^, which contain ft-quark(s) 
(much heavier than c-quark), only the 'gluon-gluon fusion' mechanism for the production 
is accurate enough, therefore, here only the 'gluon-gluon fusion' mechanism is available. 
Moreover, since the polarization of the double-heavy baryons is also strongly effected by 
hadronization of the double-heavy diquark produced via the mechanisms considered here, 
so in the present generator only the un-polarized production for the baryons are available. 

Typical running time: It depends on which option one chooses to match PYTHIA when 
generating the events and also on which mechanism is chosen for generating the events. 
Typically, for the most complicated case via gluon-gluon mechanism to generate the mixed 
events via the intermediate diquark in (cc) p.Sijg and (cc)[ 1 S'o]6 states, then on a 1.8 GHz 
Intel P4-processor PC-machine, if taking IDWTUP=1 for PYTHIA option (the meaning 
will be explained later on), it takes about 20 hours to generate 1000 events, whereas, if 
IDWTUP=3 (the meaning will be explained later on), it takes only about 40 minutes to 
generate 10 6 events. 
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Keywords : Event generator; Hadronic production; Double-heavy baryon (S cc , E bc , S&&). 



I. INTRODUCTION 

SELEX Collaboration has reported their observation of the doubly charmed baryon 

0,3 

and their measured decay width and production rate are much larger than those 
predicted theoretically, e.g. Refs. 0, 4]. Nevertheless, it started a new stage of the double- 
heavy baryon study. The observations of the double- heavy baryons Hqq/, Q,Q' = b,c 2 should 
be done in all possible experimental equipments, e.g. the detectors at the hadronic colliders 
(TEVATRON and LHC) and fixed target ones as well etc. Especially, for the lightest one 
S cc so great discrepancy between the experimental data and the theoretical predications 
should be explained. Theoretically in order to explain the discrepacy, more configurations 
of the double-heavy diquark inside the baryons and more mechanisms, in addition to the 
gluon-gluon fusion mechanism, have been taken into account (5], Q] . Considering the usages of 
experimental feasibility studies and the fact that the efficiency to generate the double-heavy 
baryon Sqq/ with PYTHIA [7] directly is too low due to the reason as that in generating 
B c meson, here in a PYTHIA-compatible format we write an effective generator GENXICC 
for the hadronic production of double-heavy baryons involving all of the formulae appearing 
in Refs. and extended to the production for the double-heavy baryons E bc and E bb 
(heavier than S cc ). Especially, the double-heavy baryon S bc has quite large branching ratio 
for decaying into J/ip inclusively, it may contribute serious background in observing B c 
experimentally. Since the format, GENXICC can be easily implemented into PYTHIA 
to simulate the full events for various experimental environments such as in the case of 
BCVEGPY for the B c meson hadroproduction , 111) . Though the double- heavy baryons 
E bc and E bb (heavier than E cc ) in the generator are involved, the cross-sections of E bc and E bb 
are much smaller than that of S cc and they are more difficult to be observed experimentally 
than S cc . To meet the most needs in cases and to make the generator GENXICC as compact 
as possible, here we only involve the most important (great) gluon-gluon fusion mechanism 



2 Throughout the paper, ^qq< denotes the baryons of (QQ'q), q = u, d a light quark, and the isospin- 
breaking is ignored. For instance, S cc denotes either S+ or S+ + etc. 
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12| . For 



for the Sfc c and E&6 production. 

The formulation for the production in Refs. is under NRQCD framework [ 
instance, the S cc -baryon production is realized via production of a binding diquark of heavy 
flavors (a diquark core) either (ccJpSijg (in configuration of S*-wave and in color anti-triplet 
3) or (cc)[ 1 S'o]6 (in configuration of S"-wave and in color sextuplet 6), and then the relevant 
double-heavy baryon is formed by means of the diquark as a core to absorb a proper light 
quark from 'vacuum' non-perturbatively. Since there is no precisely theoretical estimate 
on the relevant non-perturbative matrix element which describes the fact of combining two 
heavy quarks into the diquark and absorbing a proper quark from 'vacuum' to form the 
baryon non-perturbatively, thus as in Refs.[3j, U, |5|, [6|, the matrix element here is approx- 
imately taken as a 'decay constant' of the heavy diquark core and the probability of the 
diquark to form the baryon is assumed to be one. Namely to study the hadronic produc- 
tion of Hqq/ is equivalent to study the hadronic production of (QQ')-diquark. It may be 
understood as that the fragmentation function D(z) of a heavy diquark into a baryon has a 
very sharp peak near z — 1 |4j, and the momentum of the final baryon may be considered 
roughly equal to the momentum of the produced diquark. 

In the generator, for S cc production we include the gluon-gluon fusion mechanism that 
via the subprocess g + g — > H cc + c + c and the extrinsic charm mechanisms those via 
the subprocess g + c — > S cc + c (gluon-charm collision) or the subprocess c + c — > S cc + g 
(charm-charm fusion), furthermore, two configurations of the binding diquark, (cc) [ 3 5i]g and 
(cc)[ 1 S'o]6, are considered as in Refs.0, [a]. In the mentioned mechanisms, the heavy quark 
components in the initial hadrons are those created from 'gluon splitting' perturbatively, so 
we call them as 'extrinsic' ones for convenience. Because it contains certain non-perturbative 
nature j^], the extrinsic charm fusion via the subprocess c + c — > S cc is not included. There 
are also non-perturbative intrinsic charm components in the initial hadrons [8 , which may 
give sizable contributions when the CM. energy of the collision is not too great 6j. Since the 
contributions from the 'intrinsic charm' mechanism are model-dependent, i.e., they depend 
on the non-perturbative intrinsic charm distribution function, so in the generator we do not 
include the 'intrinsic' mechanism. 

The technology, which we will adopt here, for simulating the production by taking the 
gluon-gluon fusion mechanism and the extrinsic charm mechanism into account is that in 
Ref.Qj. Namely these mechanisms are consistently dealt with under the general-mass flavor- 
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number (GM-VFN) scheme jl^, Since here we will take the way to deal with the most 



complicated mechanism of gluon-gluon 



id, 
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has been described quite well in Refs.jsl 
it here (the interesting reader may consult Ref. 
at due places. 



usion^as that of the generator BCVEGPY, which 
so to shorten the paper, we will not repeat 
|), but we will highlight only the differences 



Due to the fact that 6-quark is much heavier than the c-quark, to produce a bb pair 
is more difficult than to produce a cc pair. Therefore, the production of E&6 and Eb c , in 
contrary to S cc , via extrinsic charm and/or bottom mechanisms is much less important 
in comparison to that via the gluon-gluon fusion mechanism. So considering the accuracy 
needed in the most usages for simulating the double-heavy baryons E&& and E^ c , we think that 
it is enough to consider the gluon-gluon fusion mechanism but not the extrinsic charm and/or 
bottom mechanisms, therefore, in the generator GENXICC, only the gluon-gluon fusion 
mechanism is involved. The GENXICC package is written in the format of PYTHIA^], so 
the generator could be easily implemented into PYTHIA. When the generator GENXICC is 
implemented in PYTHIA, one can generate complete events of the production conveniently 
as long as the created quark(s) and gluon in the hard subprocess are connected to the 
relevant subprogram for their fragmentation in PYTHIA. For user's convenience, we also set 
up a switch in GENXICC to decide whether to use the VEGAS package [3] or not, that 
depends on whether user would like to increase the Monte Carlo (MC) simulation efficiency 
for high dimensional phase-space integration with VEGAS functions or not. If using the 
VEGAS package, one may obtain the sampling importance function that is necessary for 
the importance sampling MC method. 

The paper is organized as follows. In Sec. II, we outline the generator GENXICC, which 
includes an explanation of the generator structure, the use of generator and some checks of 
the generator. The final section is reserved for a summary. 



II. THE GENERATOR: GENXICC 



Program GENXICC is a specific generator for hadronic production of a baryon S cc , S;, c 
and S&b in a Fortran package form. Since the hadronic production of E&& and Sf, c are similar 
to the case of S cc , in the paper we take the hadronic production of H cc as the main clue to 
explain the generator GENXICC, and then share some lights on the production of Sf, c and 
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FIG. 1: The schematic Feynman diagrams for the hadroproduction of H cc from the gluon-gluon 
mechanism, where the dashed box stands for the hard interaction kernel. k\ and &2 are two 
momenta for the initial gluons, q c i and q C 4 are the momenta for the two outgoing c, P is the 
momentum of H cc . The (cc)-diquark pair is either in (cc)g [ 3 Si] or in (cc^^So] respectively. 

E bb in suitable places. 

To produce S cc by hadronic collisions, the comparatively important mechanisms have: 
the gluon-gluon fusion one via the subprocess g + g — > S cc + c + c and the extrinsic charm 
one via the subprocess of the gluon-charm collision g + c — > H cc + c or the charm-charm 
collision with a hard gluon emission c + c — > S cc + (7. The schematic Feynman diagrams 
for the gluon-gluon fusion mechanism are shown in Fig.([T]). Fig.([T]) shows the two ways to 
form the (cc)-diquark pair from the outgoing valence quarks, and each way simply contains 
36 Feynman diagrams that are similar to the case of hadronic production of B c (all the 
diagrams can be found in Ref. 9|, but one need to change all the b quark line there to the c 
quark line) . For each of the sub-processes g + c — > S cc + c and c + c —>■ H cc +g of the extrinsic 
charm mechanism, there are ten Feynman diagrams, and the typical Feynman diagrams of 
them are shown in Fig.(T5]). Moreover, for every mechanism the intermediate diquark may 
be in the configurations either (cc) [ 3 Si\^ or (cc)[ 1 S , ]6- In the generator, the program is 
written as that the events are generated by specific mechanism and configuration of the 
double-heavy diquark separately, whereas if it is needed to simulate the true situation, the 
events generated by these distinguishable mechanisms and diquark configurations can be 
'mixed' by taking proper options provided by PYTHIA. More exp 
mixing events with PYTHIA can be found in PYTHIA mannual 

Since GENXICC is written in the format as that of PYTHIA (including common block 
variables), so all of the mentioned mechanisms and double- heavy diquark configurations 
can be easily implemented into PYTHIA as an external process according to user's specific 
wish. Thus all the functions in PYTHIA can be applied to GENXICC conveniently. We 



anation about generating 

3- 
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(f) (g) (h) (i) (J) 

FIG. 2: Typical Feynman diagrams for the sub-processes induced by 'extrinsic' charm. The upper 
five of them are those for g(jp\) + c(j>2) — > r. C c{pz) + c(p4) and the lower five are for c(pi) + c(p2) — > 
^cc(p3)+g(p4) respectively. The (cc)-diquark pair is either in (cc)g [ 3 Si] or in (cc)e[ 1 So] respectively. 

should note here that one has to set the PYTHIA parameter mstp(61) = when an extrinsic 
massive charm mechanism is considered, that means the showers for the initial state particles 
especially for that of the massive c-quark are switched off. It is because that the present 
version of PYTHIA can deal with the initial state showers only for massless particles 3 . 

GENXICC may be more conveniently used to simulate the events according to various 
experimental situations with a GNU C compiler under the LINUX environment at a PC 
than anyone under the WINDOWS operating system. Thus we suggest that one had better 
to work with GENXICC under a LINUX environment with FORTRAN 90 (or FORTRAN 
77) and GNU C compiler being installed. The command make in GNU is used to arrange 
the whole execution of the generator. When running the GNU command make, firstly all 
the C/C++ codes in the configuration file makefile, which is in the main directory, are 
compiled sequentially, and all of the FORTRAN source files of the program are followed to 
be also compiled, then all of the compiled files are linked according to these C/C++ codes 
automatically, finally an executable file with the default name run will be built. By running 
the built executable file run, one can generate the double-heavy baryon S cc (S fec , E bb ) events 
accordingly. Under this way, the time for compiling the Fortran source files is saved a lot. 

3 In PYTHIA, one way to solve the problem is introduced by transforming the initial massive particle to 
an equivalent massless one 0, but it fails for the present case. 
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This point will be explained more in the following subsection. We note here that all the 
sources files in our package are written in FORTRAN only i.e. its C/C++ version is not 
available so far, while as described above, alternatively one may use the GNU C compiler 
make to run the generator with FORTRAN source files to meet all kinds of needs instead. 
Furthermore, a simple script, named as do, which does all the necessary jobs for generating 
the events, e.g. calling the command make to generate the executable file run and running 
the builded executable file run and etc. is put in the main directory of the program. 

We also note that the polarization of the production for the double-heavy baryons has 
not been programmed. In fact, to estimate the polarization of the produced double-heavy 
baryons is much more complicated problem than that of unpolarized one. At this moment, 
in the program we, also as the authors in literature, treat the production in such a way just 
according to NRQCD factorization: if taking S cc as an example, firstly a proper double-heavy 
diquark (cc) in 1 S'o or 3 Si configuration is produced, that can be computed by pQCD, then 
the diquark to form the baryons in terms of non-perturbative QCD, that is not calculable 
and is dictated by the so-called 'matrix element'. Furthermore as that in literature to 
deduce the input parameters as possible, the polarization of the baryon is summed up in the 
program for the production via the diquark and spin-symmetry for the matrix element is 
assumed, that the 3 Si diquark configuration to produce a double-heavy baryon in spin-half 
or spin three halves by capture a light quark from vacuum is 'averaged' respect to spin of 
the binding diquark. Therefore, there is a weight factor, |, for producing the double-heavy 
S* c in spin three halves and a weight factor, |, for producing the double- heavy E cc in spin 
a half with the assumption. Since the excited state S* c (spin in three halves) will decay 
into the ground state (spin in one half) with an almost 100% possibility via strong and/or 
electromagnetic interaction, so considering 'comparatively poor' abilities of experiments in 
the foreseen future and for simplifying the problem, we consider the production of S* c and 
S cc together, so that without distinguishing the two states the weight factors of them should 
be added (to become 1). In this sense later on, if there is no special emphasis, the production 
of E cc means the production of S cc and S* c as well always so the matrix elements hf and 
hf mean the hadronization of the diquarks to the baryons H cc and S* c as well. For the 
production of and here the treatments are the same. Even so, for the production 
concerned, we still need to know the values of the matrix elements hf, hf, h b 3 b , hf, h b 3 c : h lbc , h\ c 
and hf c due to the various intermediate double-heavy diquarks. If someone is interested in 



9 



genxicc/ 



main directory 

phase/ generating allowed phase space 

generate/ generating events with other modules 

system/ recording running information 

pybook/ using PYBOOK to record data 

ggsub/ calculating gg-subprocess square of amplitude 

gcsub/ calculating gc-subprocess square of amplitude 

ccsub/ calculating cc-subprocess square of amplitude 

pythia/ PYTHIA 6.409 and its related subroutines 

data/ all the obtained data files are put here 



FIG. 3: The schematic structure of GENXICC. 

polarization of the 3 S\ configuration specially, only for the gluon-gluon fusion mechanism 
he can realize his interest by making proper changes in GENXICC, i.e., to make suitable 
changes in the subroutine bundhelicity (in the file s_bound.f of the subdirectory ggsub) 
accordingly, he can generate the diquark with a definite polarization as he likes. More 
explicitly, for GENXICC in default status, summation over the three components of the 3 Si 
diquark's polarization is programmed in the subroutine bundhelicity, therefore, the one 
can achieve the polarized result by simply keeping the terms corresponding to the component 
of the diquark's polarization of his interest, and commenting out the remaining terms related 
to his not interest components of the diquark's polarization in bundhelicity. 



A. Structure of the generator 

The schematic structure of the program is shown in FIGJ3J It shows that there are totally 
nine subdirectories in the main directory genxicc. The subdirectory data is used to store 
all the obtained data files, while the other eight subdirectories contain necessary files to 
complete specific tasks accordingly for generating events. Namely the generator GENXICC 
has been divided into eight modules with less cross communication among the modules, and 
each module is applied to complete assigned task(s). In the main directory genxicc, there 
are three Fortran source files: parameter. F, run.F and xicc.F. The file xicc.F is the main 
program of the generator, which arranges the running of all the Fortran codes in the main 
directory and in subdirectories in proper order. The file run.F is to define a few of global 
variables (all of which are named as the variables in capital, e.g. ENERGYOFLHC (the 
collision energy for LHC) and NUMOFEVENTS (the number of events to be generated) 
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etc. The file parameter. F is to set the initial values for all the other parameters, such as the 
mass of c-quark etc.. The eight modules of GENXICC include: 



The module generate: it is the key module of the program and it contains the files 
for generating E cc , Ef, c and events respectively with the help of the modules phase, 
ggsub, gcsub, ccsub and pythia. The functions of the module generate are to set 
the initialize condition for event simulation: to establish the connection with PYTHIA 
and a connection to a specific version of the parton distribution functions (PDFs) 
that have not been included in PYTHIA yet but included in t he p ackage for conve- 



nience (i.e. the three types of PDFs: CTEQ6HQ [15], GRV98L ll6| and MRST2001L 



17]] ) ; to calculate the integrand (the squared amplitude with necessary Jacobian for 
the phase space integration) with the help of either the module ggsub for gluon- 
gluon fusion mechanism or the module gcsub for gluon extrinsic-charm mechanism 
or the module ccsub for gluon extrinsic-bottom mechanism; to do the phase space 
integration with the help of the module phase. The module generate contains five 
Fortran source files: evntinit.F, outerpdf.F, totfun.F, sepfun.F and mixfun.F. Particu- 
larly, for simulating Sqq' production with a given mechanism, sepfun.F is to calculate 
the integrand (the squared absolute amplitude with Jacobi determinant of the phase 
integration) for a certain (QQ')~ciiquark production channel as specified by the pa- 
rameter ixiccstate (the meaning of which will be explained in the next subsection); 
mixfun.F is to calculate the integrand for the summed results, e.g. to calculate the 
summed results of the (QQ')-diquark in possible spin-color configurations for a certain 
mechanism as specified by the parameter imixtype (the meaning of which will be 
explained in the next subsection). 

The module phase: it contains the files for generating the allowed phase-space point 
and for generating an importance sampling function with the help of the VEGAS 
program fla ]. There are four source files: phase_gen.F, phase_point.F, charmsub.F 
and vegas.F. The file phase_point.F is to generate the phase-space points allowed 
kinematically with the help of the file phase_gen.F; the file phase_gen.F contains the 
RAMBO subroutine |l^ |. which can be used to generate the allowed phase space and 
to generate the momenta of the particles in final state; the file charmsub.F is only 
used in the case of S cc production for calculating the subtraction terms of the charm 
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PDF, that is necessary to avoid the 'double counting' when summing up the results 
of the two mechanisms: the gluon-gluon fusion and the gluon-charm collision [5]; the 
file vegas.F contains the VEGAS program, which is used to generate the importance 
sampling function for MC simulation. 

The module ggsub: it contains the files for calculating the squared amplitudes of E cc , 
H&& and Eb c production through the gluon-gluon fusion subprocess. As for E cc or Ebb 
simulation, when computing the squared amplitudes, the intermediate diquark is in 
either one of the two spin-color configurations (QQ)[ 3 Si\^ and (QQ)[ 1 S ]q (Q = b,c) 
respectively. While for Eb c , the intermediate diquark can be in either one of the four 
spin-color configurations (bc)[ 3 Si\^, (bc)[ 3 Si]Q, (6c)! 1 So] 3 and (fee)! 1 So] 6- There are five 
Fortran source files: s_bound.F, s_common.F, s_foursets.F, s_free.F and s_samp.F. The 



files are taken from BCVEGPY 11], but necessary changes are made so as to make 
them suitable for the present purpose (i.e., to compute the E cc , E& c and Ebb production). 
Here due to the replacement of an antiquark by a quark, the changes are included in 
color factors, the connection between the free quark part amplitudes and the bound 
state part amplitudes and etc.. More differences between the production of E cc (E bc 
and Ebb) and the production of B c through the gluon-gluon fusion subprocess can be 
found in Ref.js]. 

The module gesub: it contains files for calculating the amplitudes squared for E cc 
production through the gluon-charm collision (the gluon extrinsic charm mechanism). 
When computing the amplitudes squared, and the intermediate (cc)-diquark configura- 
tion either (cc)[ 3 Si]3 or (cc)[ 1 So]6 is considered respectively. There are two source files: 
gcampl.F and gcamp2.F. The file gcampl.F is to calculate the amplitude squared for 
the configuration (cc) pSjg and the file gcamp2.F is to calculate the amplitude squared 
for the configuration (cc)[ 1 S ]6- 

The module ccsub: it contains the files for calculating the amplitudes squared for 
producing E cc through the charm-charm collision (via the subprocess c + c — > E cc + g) . 
When computing the amplitudes squared, the intermediate (cc)-diquark configuration 
either (cc)[ 3 Si]g or (cc)[ 1 So]6 is considered respectively. There are two source files: 
ccampl.F and ccamp2.F. The file ccampl.F is to calculate the squared amplitude for 
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the configuration (cc) [ 3 Si]s and the file ccamp2.F is to calculate the amplitude squared 
for the configuration (cc)[ 1 S , ]6- 

• The module pythia: it contains the files for PYTHIA running. There are four source 
files: py6409.F (a nickname for PYTHIA 6.409), upinit.F, upevnt.F and colorflow.F. 
Here upinit.F and upevnt.F are two files to initialize an external process defined 
by user for PYTHIA. While colorflow.F sets the color flow information for the con- 
cerned processes. For convenience, the manual (pythia_manual.tar.gz) and some notes 
(pythia6409. update) for the current PYTHIA version 6.409 are included. 

• The module pybook: it contains the files for initializing the subroutine PYBOOK 
of PYTHIA to record the events. The subroutine PYBOOK can be easily adapted 
to achieve the differential distributions of a certain particle in final state. Note here 
that one can record all the necessary information of the events at one run by properly 
using of PYBOOK. Some typical ways to record the events with PYBOOK have been 
programmed in the generator GENXICC. Furthermore, the user may conveniently 
switch off this module in xicc.F in main directory so as to use his/her own way to 
record the data. There are five source files: pybookinit.F, uphistrange.F, uppydump.F, 
uppyfact.F and uppyfill.F. Here, the four files: pybookinit.F, uppydump.F, uppyfact.F 
and uppyfill.F are to call the PYTHIA subroutines PYBOOK, PYDUMP, PYFACT, 
PYFILL respectively; uphistrange.F is to set the ranges for the histograms. 

• The module system: it contains the files to open or to close the relevant record files 
and to print out certain running messages at the intermediate steps according to one's 
wish. They tell the user at which step the program is running. There are seven source 
files: upopenfile.F, uplogo.F, vegaslogo.F, updatafile.F, upclosegradefile.F, uperror.F 
and upclosepyfile.F. Here upopenfile.F and upclosegradefile.F are used to open or to 
close the files for recording the sampling importance function; updatafile.F and up- 
closepyfile.F are used to open or to close the data files used for PYBOOK; uperror.F is 
to print some possible warning messages during program running; uplogo.F is to print 
the logo of the generator to the screen or to the data file with suffix '.cs'; vegaslogo.F 
is to print the running information during VEGAS running to the screen or to the 
data file with suffix '.cs'. 
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All the files for recording the running information are put in the subdirectory data. 
To distinguish them, all the grade files (sampling importance function) for recording the 
obtained data are ended with the suffix '.grid', all the intermediate files, which record the 
used parameter values and the VEGAS running information, are ended with the suffix '.cs' 
and all the files which record the differential distributions, e.g. the transverse momentum 
and rapidity of a specific particle in final state, such ended with the suffix '.dat'. 

Each module is equipped with its own makefile that will be used to make a library of 
the same name, e.g. the makefile in the subdirectory generate will be used by the GNU 
command make to generate a library named as generate. a, which is located in the main 
directory. These sub-makefiles are orchestrated by a master makefile in the main directory. 
Under this way, the time for compiling Fortran source files can be saved a lot, because once 
a source file is compiled, one does not need to re-compile it again, unless some changes are 
made. Libraries required for the main program are listed in the LIBS variable of the master 
makefile and built automatically by invoking the sub-makefiles: 

LIBS = generate.a phase. a system. a pybook.a pythia.a ggsub.a gcsub.a ccsub.a 

Under the way based on makefile, with the command make an executable file is built 
whose default name is run, then the program acquires good modularity and re-usability. 
Thus, the user can easily reform the generator to suit his/her experimental environment. A 
simple script, which is named as do and does all the necessary jobs for generating events, is 
put in the main directory. For convenience, another script, taken directly from the FormCalc 
package [2c|, is also supplied, which is named as pnuglot and may produce a high-quality 
Encapsulated PostScript (EPS) file to plot the generated data in relevant file(s) that is(are) 
in the subdirectory data. 

B. Flow charts of generator 

The overall schematic flow chart of the generator GENX1CC is shown in Fig.pj. It 
shows that GENXICC can be schematically divided into two blocks, i.e. the vegas block 
(mainly with phase) and the event block (mainly with generate). The function of the 
vegas block is to generate the sampling importance function that is necessary for importance 
sampling, and the function of the event block is to generate events with PYTH1A, e.g. all 
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FIG. 4: The overall schematic flow chart of the generator GENXICC. 
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FIG. 5: The schematic flow charts for the vegas block (on the left) and the event block (on the 
right) of the generator GENXICC. 



the mentioned production processes for the hadornic production of z<qq' can be implemented 
straightforward in PYTHIA as external process by properly programming the two PYTHIA 
subroutines UPINIT and UPEVNT 4 . The schematic flow charts for the vegas block and the 



4 The UPINIT and UPEVNT are two necessary subroutines to define external process for PYTHIA: UPINIT 
is called by PYTHIA subroutine PYINIT and is used to set the initial conditions for the user-defined 
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event block of the generator GENXICC are shown in Fig.flSj). 

When simulating the production, one may chose whether to take importance sampling 
method with VEGAS or not. If the importance sampling method is used the simulation 
efficiency may increase a lot. The VEGAS package is programmed to achieve a sampling 
importance function. One can learn from Fig.(j3J) that there are three ways to do the event 
simulation in the program: one is the trivial MC and the left two ways are to use the 
importance sampling method, i.e. one is to use the existent grades (importance sampling 
function) that have been generated in previous run(s) and have been recorded in .grid file(s) 
in the subdirectory data and the other one is to use a new grade generated by the current 
VEGAS running. When setting IVEGASOPEN=0 and IGRADE=1 for instance when the 
function has been recorded in .grid file(s) by previous run(s) already, i.e. adopting the first 
importance sampling method, one may generate S cc events, just by reading the existent 
importance sampling function. Thus, when using the package under these proper options, 
one just needs to run VEGAS once enough. 

The precision of the generated importance sampling function can be improved by properly 
adjusting the maximum iteration number, the number of calls to the integrand in each time of 
the iteration, and the number of bins (the [0, 1] region is divided into how many sub-regions) 
as well. For convenience, we define three parameters in the head file invegas.h: NVEGITMX 
(maximum number of allowed iterations), NVEGCALL (maximum total number of the times 
to call the integrand in each iteration) and NVEGBIN (the number of bins). A subtle point 
is that, in VEGAS, the default number of NVEGBIN is 50 (i.e. the region of (0,1) is 
divided into 50 pieces); by taking a bigger proper value for NVEGBIN, one can improve 
the production efficiency but need more iterations to obtain a stable result. In practice, the 
values for NVEGITMX, NVEGCALL and NVEGBIN should be carefully chosen to obtain 
the best important sampling function within the least of cpu time. In our program, we take 
the default value of the column bins to be 100 (according to our experience, it is enough); 
and the user who wants to do some very precise studies, a proper variation of NVEGBIN 
is needed. And in practice, we have found that for the most complex gluon-gluon fusion 
mechanism for the hadronic production of S cc , and Sfe c , if choosing NVEGITMX=20 

process, and UPEVNT is called by PYTHIA subroutine PYEVNT and is used to set the event conditions 
of the user-defined process. 



16 



and NVEGCALL=300000, then we can achieve the precision goal of 0.1% for the total 
cross-section. 

The flow chart for the event block is shown in the right diagram of Fig. ([5]). It is to 
generate events with the help of PYTHIA package. At the present, we adopt the PYBOOK 
subroutines that are provided by PYTHIA, such as PYFILL, PYDUMP, PYFACT and etc., 
to record all the obtained information of the events. However, users may comment out all 
the calls for PYBOOK subroutines in the source file xicc.F and then add his/her own ways 
to record the events conveniently. 



C. Use of the generator 

GENXICC can generate a huge event sample of producing S cc or E bc or E bb efficiently 
in the case of hadronic collisions. As for the hadronic production of Sno (Q = c,b), in 
f _ k of NRQCD H, t he parted ™ctt» o f I _d ZL* 1 b e 

nn 

expressed as [6|, |21[ : 

a ab ^ QQ = H(ab - (QQ)[ 3 Si]~ 3 ) • h$ Q + H(ab - (QQ^SoU) • h? Q + • • • , (1) 

where the ellipsis stands for the terms in higher orders of v. H(ab — > {QQ)[ 3 Si] 3 ) or 
H(ab —* (QQ)l 1 So]e) is the perturbative coefficient for producing a QQ pair in configuration 
of 3 Si and color 3, or in configuration of 1 5'o and color 6 respectively. The parameters 
and hf are the values of the nonperturbative 'matrix elements' in NRQCD to characterize 
the transitions of the charm pair of quarks in [ 3 Si]3 and in [""So] 6 into the baryon Hqq 
respectively. According to the argument presented in the INTRODUCTION, the hadronic 
production of Eqq is 'equivalent' to the hadronic production of (QQ)-diquark. Under such 
condition, the value of NRQCD matrix element can be related to the wavefunction for 
the color anti-triplet [ 3 Si] (QQ)-diquark state, i.e. h® Q = 5 ■ |-Rqq(0)| 2 /47t, where Rqq(0) is 
the wave function at origin for the diquark binding system (QQ) and 5 is the weight factor 
from 'polarization average assumption' in producing the baryon that is mentioned above). 
While as for Sb c , the factorization formula is different from that for Sqq due to different 
flavors being involved and, in fact, there are two physical states Ei bc and E 2bc , thus it reads: 

t ab ^s lbc = H(ab - (bc)[ 3 Si]- 3 ) ■ hf + H(ab - (6c)[%] 6 ) • h\ c + • • • 
o a ^ hc = H(ab - (bc^Sok) ■ hf + H(ab ^ (bc^Sfo) ■ h'\ c + ■■■. (2) 
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TABLE I: All the considered channels for the hadronic production of the double-heavy baryon, 
which are defined by the two parameters mgenxi and ixiccstate. Here the symbol gg-channgel 
stands for the gluon-gluon fusion channel and etc.. 





mgenxi=l (for H cc ) 


mgenxi=2 (for Hf, c ) 


mgenxi=3 (for Ebb) 


ixiccstate=l 


gg-channel, (cc)g ( 3 Si) 


gg-channel, (60)3 ( 3 Si) 


gg-channel, (66) 3 ( 3 Si) 


ixiccstate=2 


gg-channel, (cc) 6 ( 1 S , ) 


gg-channel, (6c) 6 ( 1 5 ) 


gg-channel, (66) 6 ( 1 5 ) 


ixiccstate=3 


gc-channel, (cc)^( 3 Si) 


gg-channel, (bcje^Si) 




ixiccstate=4 


gc-channel, (cc)e( 1 5'o) 


gg-channel, (6c)g( 1 5 ) 




ixiccstate=5 


cc-channel, (cc)g( 3 5i) 






ixiccstate=6 


cc-channel, (ccjei^So) 







More explicitly, four matrix elements: h\ c , h'\ c , h^ c and h'^, which characterize the transi- 
tions of the bottom-charm diquarks in pSoJe, [ 3 Si]6, [ 3 *S , i]g and pSojg into Eb c respectively 
are needed. Similarly /z| c can be related to the wavefunction of the color anti-triplet diquark 
(&c)[ 3 Si] as h b 3 c = 5 ■ |i? 6c (0)| 2 /4vr. 

Users may communicate with or give instructions to the program through the source files 
run.F and parameter. F. These two input files allow users to setup the generation parameters 
and requests. All the necessary parameters are: 

• pmc=: mass of c quark (in units GeV). It should be taken as an effective one that can 



be derived from the value of the pole mass and the MS running mass |22j. Its default 
value is taken to be 1.75 GeV. 

pmb=: mass of b quark (in units GeV). Similarly, it should be taken as an effective 
one and its default value is taken to be 5.1 GeV. 

mgenxi=: symbol to choose which double-heavy baron, i.e. S cc , and Sfe c , to be 
generated. mgenxi=l is to generate S cc events, mgenxi=2 is to generate Eb c events 
and mgenxi=3 is to generate "Ebb events. As for the production of and Sb c , only 
their dominant gluon-gluon fusion mechanism is programmed. For clarity, we show all 
the considered channels in TABJU In default, we set mgenxi=l. 
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pmxicc=: the mass of H cc , Eb c and E&& (in units GeV) respectively according to 
the value of mgenxi: pmxicc=M= cc for mgenxi=l; pmxicc=MH bc for mgenxi=2; 
pmxicc=MH bb for mgenxi=3. To ensure the gauge invariance of the hard scatter- 
ing amplitude, we must set pmxicc=pmc+pmc for S cc , pmxicc=pmb+pmc for Sb c and 
pmxicc=pmb+pmb for 

fxicc=: the value of R cc (0) or i2& c (0) or Rbb{0) (in units GeV 3 / 2 ) corresponding 
to the value of mgenxi: fxicc=i? cc (0) for mgenxi=l; fxicc=i?b c (0) for mgenxi=2; 



fxicc=i?;,fe(0) for mgenxi=3. The defau 



0.70GeV 3/2 , 0.90GeV 3/2 and 1.38GeV 3/2 js] respectively. 

cmfactor: the relation among the non-perturbative matrix elements as defined above. 
It has been pointed out in [2l| that the relevant matrix elements have the values at 
the same order of v. So for convenience, we introduce an addition parameter cmfactor 
which may vary in a possible region and further assume that ~ cmfactor • hf, 
h bb ~ cmfactor ■ h b ^ and h 1 ^ ~ h /bc ~ h' h £ ~ cmfactor ■ h b £. One may study certain 
dependence of the production on the matrix elements by setting the value of the 
parameter. Note that here we take cmfactor = 1.0 as its default value. 

ptcut=: (in units GeV) cut for S cc , E bc and E bb respectively. In default, we set 
ptcut=0.2GeV. 

etacut=: rapidity cut for S cc , E bc and E b b respectively. In default, there is no rapidity 
cut for the production. 

psetacut=: pseudo-rapidity cut for S cc , Eb c and respectively. In default, there is 
no pseudo-rapidity cut for the production. 

npdfu=: choice of the collision type of hadrons. npdfu=l for p — p collision with 
ecm = 1.96 TeV (TEVATRON), npdfu=2 for p — p collision with ecm = 14.0 TeV 
(LHC) and npdfu=3 for a fixed target experiment, but in the program we take it such 
as p — p collision with ecm = 33.56 GeV (SELEX). In default, we set npdfu=l. 

ecm=: total energy for the hadron collision (in units GeV). It is set to be ENER- 
GYOFTEVA, ENERGYOFLHC or ENERGYOFSELEX in accordance to the value 
of the parameter npdfu. 
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t values for mgenxi =1,2, 3 are taken as 



• imix=: whether to generate the mixed events. imix=0 for generating separate events 
(e.g. events from a particular mechanism and a fixed configuration of the intermediate 
diquark state) according to the value of ixiccstate; imix=l for generating mixed events 
according to the value of imixtype. As for mgenxi=l, it is to generate mixed events 
of (cc) pSiJjj and (cc) [ 1 £'o]6 of the same mechanism. As for mgenxi=2, it is to generate 
mixed events of (bc)[ 3 Si]s, (6c)[ 3 S , i] 6 , (6c)[ 1 S , o]g and (6c)[ 1 S , ]6 of the gluon-gluon fusion 
mechanism. As for mgenxi=3, it is to generate mixed events of (bb)[ 3 Si]a and (66)[ 1 5 , ]6 
of the gluon-gluon fusion mechanism. In default, we set imix=0. 

• ixiccstate=: symbol for producing the events from different mechanisms and different 
intermediate diquark states in the case of imix=0. When it is used for the production of 
H cc , the option of ixiccstate is G [1, 6] : ixiccstate=l stands for the case of gluon-gluon 
fusion with the intermediate diquark in (cc) [ 3 Si]s state; ixiccstate=2 stands for the case 
of gluon-gluon fusion with the intermediate diquark in (cc) [ 1 S'o]6 state; ixiccstate=3 
stands for the case of gluon-charm collision with the intermediate diquark in (cc)[ 3 «Si]3 
state; ixiccstate=4 stands for the case of gluon-charm collision with the intermediate 
diquark in (cc)f 1 S ] e state; ixiccstate=5 stands for the case of charm-charm collision 
with the intermediate diquark in (cc)[ 3 S , i]g state; ixiccstate=6 stands for the case of 
charm-charm collision with the intermediate diquark in (cc)[ 1 S'o]6 state. When it is 
used for the production of E bb , the option of ixiccstate is G [1,2]: ixiccstate=l stands 
for the case of gluon-gluon fusion with the intermediate diquark in (bb)[ 3 Si]s state; 
ixiccstate=2 stands for the case of gluon-gluon fusion with the intermediate diquark 
in (66)[ 1 S'o]6 state. When it is used for the production of S fec , the option of ixiccstate is 
G [1,4] : ixiccstate=l stands for the case of gluon-gluon fusion with the intermediate 
diquark in (be) [ 3 Si]s state (for Si& c ); ixiccstate=2 stands for the case of gluon-gluon 
fusion with the intermediate diquark in (be) I 1 So] q state (for E.i bc ); ixiccstate=3 stands 
for the case of gluon-gluon fusion with the intermediate diquark in (&c)[ 3 5i]6 state 
(for S 2 fe c ); ixiccstate=4 stands for the case of gluon-gluon fusion with the intermediate 
diquark in (be) pSoJg state (for H 2 b c ). And its value should be set according to one's 
needs. 

• imixtype=: the type of mixed double-heavy baryon (S cc , E bc , E bb ) events for the same 
mechanism. It is used when imix=l. In the case of S cc , imixtype G [1,3]: imixtype=l 



20 



is to generate the events for the gluon-gluon fusion mechanism with the summed up 
results of the intermediate diquark in (cc) [ 3 Si]s and (cc) pSble states; imixtype=2 is 
to generate events for the gluon-charm collision mechanism with the results of the 
intermediate diquark in (cc) [ 3 <Si]jj and (cc)[ 1 S'o]6 states summed up; imixtype=3 is 
to generate events for the charm-charm collision mechanism with the results of the 
intermediate diquark in (cc)[ 3 S , i]g and (cc)[ 1 S'o]6 states summed up. In the case of Eb c 
or S(,b, only one mixed type of event is programmed, i.e. imixtype=l is to generate 
the events for the gluon-gluon fusion mechanism with the results of the intermediate 
diquark in (bc)[ 3 Si]s, (bc)^S ]e, (bc)[ 3 Si} 6 and (be) ^So]^ states for the production of 
Sfe c summed up or those with the intermediate diquark (bb)[ 3 Si]s and (66)[ 1 S'o]6 states 
for the production of E bb summed up. And its value should be set according to one's 
needs. 

ivegasopen=: (status of VEGAS subroutine) ivegasopen=l means the VEGAS sub- 
routine is 'on'; ivegasopen=0 means the VEGAS subroutine is 'off'. In default, we set 
ivegasopen=l. 

igrade=: whether to use the grade generated by previous VEGAS running, which is 
used only in the case of ivegasopen=0. igrade=l means to use; igrade=0 means not 
to use. Thus one runs VEGAS once enough, since one can use the existed grade to 
generate events by setting ivegasopen=0 (without running VEGAS) and igrade=l. In 
default, we set igrade=0. 

iveggrade=: whether to improve the existed grade generated by the previous (earlier) 
VEGAS running in the case ivegasopen=l. igrade=l means to use; igrade=0 means 
not to use. By setting igrade=l, one can generate a more precise grade based on the 
existed grade. Such option is used in the cases when one does not satisfy with the MC 
precision of the old existed grade. In default, we set iveggrade=0. 

number=: total number of times for calling the integrand. The parameter is needed 
only when ivegasopen=l. Its value, together with the value of the following parameter 
nitmx, should be adjusted according to one's precision goal. 

nitmx=: upper limit for the number of iterations. The parameter is needed only when 
ivegasopen=l. 
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nev=: number of the events to be generated for the hadronic production of h + h — > 
S cc + • ■ •, h + h — > Ebc + ■ ■ • or h + h — > + ■ • • (h means a hadron) respectively 
according to the value of mgenxi. Its value should be set according to one's needs. 



ioutpdf=: indicates whether the three latest versions of PDFs: CTEQ6HQ [15] . 
GRV98L Q and MRST2001L [17] are used when generating the events. The three 
PDFs are offered in the program for the hadronic production. ioutpdf=l means they 
are used with further option ipdfnum, while ioutpdf=0 means they are not used so the 
inner PDFs of PYTHIA are used by the option mstp(51). In default, we set ioutpdf=l. 

mstp(51)=: choice of the proton parton-distribution set (PYTHIA parameter, see 
PYTHIA manual), e.g. mstp(51)=2 for CTEQ3M, MSTP(51)=7 for CTEQ5L, 
mstp(51)=8 for CTEQ5M, and etc.. It is used when ioutpdf=0 and in default, we 
set mstp(51)=7. 

ipdfnum=: indicates which one of the three latest version of the PDFs CTEQ6HQ 
151 ]. GRV98L [16( and MRST2001L [17J is used for the hadronic production. It comes 
into operation only under the option ioutpdf=l. ipdfnum=100 means to use GRV98L; 
ipdfnum=200 means to use MRST2001L; ipdfnum=300 means to use CTEQ6HQ. It 
should be noted that only CTEQ6HQ is consistent with the calculation technology of 
GM-VFN scheme, which is our default choice for PDF. The other two PDFs GRV98L 
and MRST2001L programmed here are only for comparison with the results in litera- 
ture easily. In default, we set ipdfnum=100. 

idwtup=: master switch dictating how the event weights and the cross-sections should 
be interpreted (PYTHIA parameter, see PYTHIA manual). When idwtup=l means, 
parton-level events have a weight at the input to PYTHIA. The events generated 
stochastically are either accepted or rejected according to the weight, so that fully 
generated events at the output have a common weight; parton-level events have a unit 
weight at the input when idwtup=3, i.e., they are always accepted 5 . In default, we 
set idwtup=3. 



For the present usages, of PYTHIA the two options idwtup=l and idwtup=3 are enough. 
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FIG. 6: The energy dependence of the integrated partonic cross-section for the production of the 
double-heavy baryons with heavy diquarks via the gluon-gluon fusion mechanism. The dotted line, 
solid line, dashed line and dash-dot line stand for the baryons with (cc)jj[ 3 Si], (bc)^[ 3 Si], (bcj^So] 
and (^7)3 [ 3 Si] respectively. Note here that the curves for H cc and E&& both are divided by a factor 
'2' for convenience in comparison with literature results. 

• mstu(lll)=: order of a s in the evaluation in PYALPS (a PYTHIA routine for calculat- 
ing a s , see PYTHIA manual); e.g., mstu(lll)=l for leading order (LO); mstu(lll)=2 
for next leading order (NLO). In default, we set mstu(lll)=l. 

• paru(lll)=: constant value of a s (see PYTHIA manual), which is used only when 
mstu(lll)=0. In default, we set paru(lll)=0.2. 

• isubonly=: whether to generate the information only of the gluon-gluon fusion sub- 
process; isubonly=0 for the full hadronic production, i.e. the structure functions are 
connected; isubonly=l for the subprocess only. In default, we set isubonly=0. 

• subenergy=: the energy (in units GeV) of the gluon-gluon fusion subprocess. It is 
needed only when isubonly=l. And its value should be set according to one's needs. 
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D. Check of the generator 



First of all, all the programs are checked by examining the gauge invariance of the am- 
plitude, i.e. the amplitude vanishes when the polarization vector of an initial/final gluon 
is substituted by the momentum vector of this gluon. Numerically, we find that the gauge 
invariance is guaranteed at the computer ability (double precision) for all these processes. 
Especially, to make sure the correctness of our program for the most complicated gluon- 
gluon fusion mechanism, we have checked whether the numerical results evaluated by the 
generator agree well with those by FDC package 0, Q]. For such purpose, we have added 
the parameter isubonly to decide whether to generate the information only of the gluon- 
gluon fusion subprocess, since it is much more easier to be compared. Furthermore, we 
have compared the numerical results for the gluon-gluon fusion mechanism with those by 
using the same input parameters as stated in the corresponding references in the literature, 
and the examples of these comparison can be explicitly found in the section III of Ref . jsj . 
In checking the program, we have recalculated all the values listed in Ref. JfJ, and we find 
a good agreement of our present results with those of Ref.jlj]. More explicitly, we obtain 
the same results for the partonic cross sections as that of Ref. J^J and we show the partonic 
cross sections for the production of baryons with heavy diquarks via the gluon-gluon fusion 
subprocess in Fig. ([6]). 

Finally, we show a test run for the hadronic production of H cc at the TEVATRON. When 
running the program, the initial parameters are shown in Fig.(JTJ). All the obtained data are 
turned into a zipped file, named as testdata.tar.gz in the main directory of the program. 
Some typical resultant curves, e.g. p t - and ^-distributions are shown in Fig.([H]). In drawing 
the distributions, several typical rapidity cut, e.g. \y cu t\ = 0.6, 0.1 and 0.2 are taken for 
^cc-Pt distributions and similarly, several typical p tcu t are taken for S cc -y distributions. 



III. SUMMARY 



A generator GENXICC has been developed and well-tested by us, which is for hadronic 
production of S cc based on three mechanisms via gluon-gluon fusion, gluon-charm collision 
and charm-charm collision. While for the hadronic production of Sb c and only the most 
important mechanism via gluon-gluon fusion is available, because we think that considering 
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initial parameters 



to generate the events for xi_cc. 
get the mixing results for xicc. 

mixed events of gg-channel of [(3 A s_1 )_cbar3+(1 A s_0)_c6] 
generate evnts 10000000; teva energy(gev) 0.196E+04 

m xicc = 3.5000 m b = 1.7500 m c = 1.7500 R cc (0)=0.70006 

Q2=square of xicc transverse mass ; alphas order=lo 

new version of outer pdfs which are from www 

pdf: cteq6hq 

using pythia idwtup= 3 

ptcut=0.000 gev 

no rapidity cut 

no pseudo-rapidity cut 

the used bin number in vegas: 100 

using vegas: number in each iter.= 500000 iter.= 20 

end of initialization 



FIG. 7: Snapshot of the initial parameters used in the test run of the generator GENXICC, which 
is for the hadronic production of H cc at the TEVATRON and is to generate the mixed events 
through gluon-gluon fusion mechanism (imix=l and imixtype=l). 




P t (GeV) 



FIG. 8: H cc -pt and H cc -y distributions for the test run of the generator GENXICC, which is for 
the hadronic production of H cc at the TEVATRON and is to generate the mixed events through 
gluon-gluon fusion mechanism (imix=l and imixtype=l). 

the characteristics and most potential usages it is enough, and the program may be simplified 
in certain degree. The generator has proper interface to PYTHIA, so the full events of the 
production can be generated by filling the standard PYTHIA event common block. Note 
that in the program we do not distinguish the state Sqq/ (in spin one half) and the state 
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Sqq, (in spin three halves) and consider their production of them together. If one would like 
to know their difference of the production of two states, roughly speaking, one can estimate 
it by the way to account the production as that if it is via 1 So heavy diquark then all (100%) 
of the diquarks will be fragmented into the state Sqq/, while if it is via 3 Si heavy diquark 
then | fraction of the diquarks will be fragmented into the state Sqq/ and | fraction will 
be fragmented into the state Hqq/. In view of the prospects for double- heavy baryon (S cc , 
S;, c and Sf>b) physics at SELEX, Tevatron and LHC, the generator offers a useful tool for 
further experimental feasibility studies of the double-heavy baryons. 
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